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DYNAMICS  OF  ANGULAR  MOVEMENTS  OF  A 
SOLID  SUPPORTING  A  ROTATING  ROTOR 
WITH  CONSIDERATION  OF  ENERGY 
DISSIPATION 

V.  A.  Grobov,  I.  I.  Kantemir 

This  article  analyzes  the  motion  of  a  system  consisting  of  an  axi- 
symmetric  solid  B  with  mass  M,  rotor  B'  with  mass  M' ,  and  two  passive 
dampers  b  and  b’,  each  of  which  consists  of  a  spring  and  mass  placed 
into  a  tube  filled  with  a  viscous  fluid.  The  schematic  of  the  system 
is  shown  in  the  figure.  Body  B  is  used  as  the  principal  body,  i.e., 
a  body  relative  to  which  the  movement  of  all  the  components  of  the  sys¬ 
tem  is  being  analyzed.  Xj,  X£,  X3  is  the  coordinate  system  coupled  with 
body  B,  the  origin  of  which  is  at  the  center  of  masses  of  body  B  and 
one  of  its  axes  -  Xg  -  is  oriented  along  the  symmetry  axis  of  the  body. 
Consequently,  axes  Xp  X2,  and  Xg  are  the  principal  central  axes  of 
inertia  of  body  B.  The  origin  of  the  X^X^Xg  system  coupled  with  rotor 
B'  is  at  the  center  of  masses  B* ’ ;  axis  X3  coincides  with  axis  X3.  The 
distance  between  B* '  and  B*  equals  /,  while  S>  is  the  angle  between  axis 
Xj  and  the  line,  which  is  parallel  to  axis  Xj  and  passes  through  point 
B*\ 

The  damping  mechanism  consists  of  mass  m,  which  moves  in  the  tube 
filled  with  viscous  fluid.  We  designate  the  coefficients  of  the  spring’s 
rigidity  and  viscous  damping  of  the  fluid  by  k  and  c.  The  tube  is  at¬ 
tached  to  body  B  in  parallel  to  axis  X3  in  such  a  way  that  the  spring 
is  not  stretched  when  particle  m  is  at  axis  X^  at  distance  a  from  point 
B*.  The  displacement  of  the  particle  relative  to  axis  Xi  is  determined 


by  value  z.  Three  additional  particles,  the  mass  of  each  of  which  is 
equal  to  the  total  mass  of  damper  mb,  are  coupled  rigidly  with  body  B 
in  such  a  way  that  as  to  have  the  system  (body,  damper,  and  three  par¬ 
ticles)  symmetrical  relative  to  axis  X3  when  the  spring  is  not  stretch¬ 
ed.  A  similar  damper  is  attached  also  to  rotor  B'  [H]. 


Fig.  1 


Let  the  <ol(  a,.  «•»  be  the  angular  velocities  of  rotation  of  body  B 
around  axes  Xi,  X2,  and  X3  in  an  inertial  frame  of  reference.  We  as¬ 
sume  that  the  angular  rotation  speed  of  the  rotor  around  axis  X3  equals 
<*>j  +  a,  where  0  =  const. 

To  .formulate  the  equations  of  motion,  first  we  will  work  out  the 
expressions  for  the  kinetic  and  potential  energy  of  the  system. 

The  kinetic  energy  of  the  system  equals 

r  “  -y-  (M  —  A*rC*  —  2AUW  +  mfi*  +  m**'s  +  2z'lmb)  (w?  +  a4)  -f 

+  C«4  +  2J3oo),  +  /jo*  +  m*  (i*  —  2a(o}a>fi  —  2 au),2)  q  ^ 

+  m*J  *'*  —  2 loV  (<o*  +  0)  (<o,  cos  $  +  <o,  sin  $)  -f  2fl'*'<i>,  sin  $  — 

—  2oV  ©,  cos  $]  —  MtC*;, 

whe  re 

A  -  U  +  /i  +  2 m*a»  +  2mio'J  +  (M'  +  4 rrib)(\  —  v)/2, 

C  «=  /*  +  /»  +  4m*a*  -f  4/nia'*, 

^3  *  /» +  4mifl'2, 

v  -  (Af '  +  4mb)/(M  +  M'  +  4m*  +  4m;), 

AfT  =  M  +  4m*  +  Af '  -f  4mi,  ( 2 > 

*'--5;.  t-^+pV, 


Ip  and  l|,  1^  are  the  principal  enertia  moments  of  body  B  and  rotor 
B* ,  respectively. 

The  potential  energy  of  the  system  has  the  form 

n  +  (3) 

Taking  into  account  expressions  (1)  and  (3),  we  write  the  equations 
of  motion  in  the  Lagrange  form: 

—  —  C)®g®,  +  J'a  o®,  +  +  Mr?  (a,  —  w^>,)  + 

+  «*!— 2  (;  +  /,)  [rftij  —  o),©,)  +  2®,]  +  2|— 2j®1  +  22®1  + 

+  2K  —  ®*“»)  —  a (a,  +  u, to,)] j  +  mb  [—  2 (5  —  /,)  \z‘ (a,  —  a,a3)  + 

+  *’®d  +  **  |—  2£®x  +  2z'®,  +  2'  (ciij  —  a^o,)  — 

—  a'cos^(«D,  +  «1oil)]  +  a'siny){'2  +  2'lK  +  o)t— c4])l  =0.  (4a) 

Ai,_(C— X)®,®,  —  J'sai^  +  2  iW^fc®,  +  +  a,®,)  + 

+  mb  { — 2(t  +  /»)lz  (®*  +  Mi®,)  +  2®,]  +  2 1 — 2j®,  +  22(1),  + 

+  *(“*  +  —  a (2  +  2 (®s  —  ©?))}  +  nib  [—  2 (£  —  /,) \z'  (a,  +  a,®,)  + 

+  z'«J  +  z'  [ —  2£<i),  22'®,  -j-  2'  (a,  -j-  a,©,)]  — 

—  a'cos^{2'+2'I(®,  +  o?  — ®?1}— a'2'sint(®»— ®i®*>] ->=0,  (46) 

(22®,  -f  2(®|  —  ©»©d] — mia'  sin  t  [2z'  •«  +  2'  (©,  +  ®,a,))  — 

—  niba'  cos  $ 122'®,  +  2'  ,w,  —  b*©,))  »  0.  (4b) 

m,(I  —  p)’z— m'bpz' — «vb  (®* — a,®,)  _  m*  (®f  +  a?)  (2  (1  —  p)  — 

—  /,  —  P'z'J  +C2  +  kz  -  0,  (5a) 

—  «*P'z  +  mi(l— p')?  +  mia'  (sin  $  I®,  +  «,  (a,  +  2o)l  — 

— eos^K  —  ®,(®,  +  2o)]}  —  mi(®?  +  *4)j2'(l  —  p')  +  /1  —  pz)  + 

+  c'i'  +  *V-0.  (56) 

Here  we  use  the  designations 

,  /(Af  +  4 ,  *<*»'+«•£) 

**“  *T  •  *»“  an - •  (6) 

We  introduce  the  following  quantity: 

,  (7) 

which  represents  an  angular  velocity  of  free  precession  of  a  system  of 
bodies.  Then  equations  (4a)— (4c)  can  be  written  as 

•S  —  A®,  -  |i F,(®i,  ®», a®  •»  2, 2,  s',  z',  z>), 

«S-  l^s  («u«W  «•*  »».«•»«.*)  *\  z'X 


where  p  is  a  small  parameter: 


and  F^,  F2,  F^  -  nonlinear  functions  in  the  form 

Fi  -  ~x  4 *»,)  —  2sa h  —  *•(«•,  —  ov*,'  + 

4  « («*  4  “H®*)}  — pr  l2ii  t*'  (°h  ~  “*®i)  4  4  ^i'a,  4 

+  *'*  (®i  —  ®*®»)  —  a'f  cos  *  («,  +  «,«»,)  +  a'  sin  \p  [?  4  ^  ((to,  4  of — 4))}, 

F*  =  — — r  (2/j2(®»  +  ®i«*)  4  «“*]  —  2zZ(Om  —  2*  (<ig  4  ffl.ttk)  4 

fli  a 

4«[*4*(®3  — «?)]) - !r(2/if2'(<a»  +  ®i«»)  +  i'o>«l+  2^?^+  ( 10) 

Q 

+  *'!(«1  +  — a’  cos  if  [*'  +  ^  ((<0,  +  ff)»  —  G)?)l  —  a v  sin  if  (to,—  , 

F*  -  C^'|'^r'(2a*»  +  +  (2*'®t  +  *'»«)  *in  *  4  (2^®1 4  «* 

( He  re  m,  mm  m  h  m»  m  m'.) 

After  introducing  the  designations 

-S;-W.  •£-*»  (ID 

-^7-  =*  <D*  —  »?  4  u»|  —  v^, 

we  write  equations  (5a)  and  (5b)  in  the  form 

i  4  2V  +  v?/  -  n^«(®*,  a  ?,  *)  4  « («* —  ®i®»). 

?42V'4-&'  -|iF,(®,,i',  'i',z')4  a'  |cos$  [®»  — ®i(®»  4  2®)1  —  ( 12) 

—  sin  f  4  «*(•»  4  2o)l). 

whe  re 

F4-M4V>4®,'(*-g.  (13) 

F»  -  h  <*  +  *')  +  «*»’(*'  +  V* 

Thus,  as  a  result  of  transformations,  system  (Ma)-(5b)  is  reduced 
to  the  system  of  equations  (8)  and  (12).  Resonance  cases  are  possible 
in  this  system,  when  one  of  the  intrinsic  frequencies  of  the  system 
coincides  with  the  frequency  of  precession.  Let  us  consider  some  of 
them,  relying  on  the  results  of  an  analysis  conducted  on  the  frequency 
spectrum,  which  are  given  in  work  [4],  where  it  is  shown  that  the  fol¬ 
lowing  resonance  relations  are  possible  in  the  system: 

Thus,  for  example,  if  the  following  relation  is  fulfilled  between 
the  parameters: 


•T“  1  + 


then  there  is  resonance  of  the  type  A  —  o. 

Let  us  stop  to  consider  one  of  these  cases.  Let  us  assume  that 
A  -  o  +  |iA.  Then  system  (8)  can  be  written  as  . 

«i  +  cHo,  fi  (Fl  —  Aw,), 

—  ow, +  (14) 

•-R/V 

Considering  the  case  of  fast  rotation  of  a  body  around  the  symmetry 
axis  X3,  i.e.,  assuming  that  it  is  possible  to  consider  «oi,o§ 

and  *»,«,  to  be  the  value  of  the  order  of  smallness  p. 

In  order  to  solve  the  system  of  equations  (14)  and  (12),  which  de¬ 
scribe  the  motion  of  this  system,  we  will  use  the  method  of  perturba¬ 
tions  [1—3].  For  this,  we  will  first  examine  an  unperturbed  system 
(when  p=0) 


•»  +  OW,  -  0, 
•H  —  «*j  “0, 

w,-0 


The  solution  of  (15)  is  found  in  the  form 


’«#,  —  <Dcos6, 
w,  =  ®  sin  6, 
,w.  «=  Q, 


whe  re  0  ■=  at  +  d. 


Using  the  solutions  of  (16)  as  the  generating  ones,  we  transform 
(14)  in  accordance  with  the  ideas  of  N.  N.  Bogolyubov's  method  of  per¬ 
turbations  [1,  2]  to  the  new  variables  Q>  and  0: 

-^(FjCose  +  Ft sin  0), 

-y-“-5-(^*cos0  — Fjsin8)  +  fiA,  ^ 17  ^ 

When  u=0,  taking  into  account  (16),  system  (12)  is  written  as 


fir  +  2h^t  +  k?*  —  fl>a(o  —  Q)cos0, 

\i*  +  2h^’  4.1^1'  —  —  ®a'(o  +  Q)cos6. 


We  introduce  the  following  new  variables  into  (18): 


5 


•-“(■Sp- +*k) 


— 0)cose, 


1*  \,3r'  +  *W«)  -  —  ®fl'(o  -f  Q)cos  #. 

where 

The  solution  for  (19)  we  seek  in  the  form 


whe  re 


Jjr,  -  A\  eat  Q,  +  //,«**  cos  (8 — pj, 
vi  ■  4w*Q|  ■+■  H/P*  cos  6, 


Q,  =  V  +  Q, -=  V  +  4-V. 

u  ®fl  (O  —  2)  .  2M 

H‘-yw=^W- 


•aMc-t-OJ 


Let  us  examine  the  perturbed  system 


e_%<  +  til/i)  =  |iFt  +  ®a  (o  —  co,)  cos  9, 

*~M  ('^L  +  bh/tj  *  »  —  <®a'  (o  +  co,)  cos  6. 


In  this  case  we  note  that  •  <D,  8, and  Hj  and  H2  will  be  variable  in 
perturbed  motion.  It  is  also  easy  to  see  that 


where 


40 

4  — 


Qi  =  — ■,  —  (Fx  cos  8  +  F,  sin  8), 

Tx 

Qt  =  — (Fx  cos 0  -f  F,  sin  8), 
'I 

R  =  A 


-The  solution  for  (21)  we  seek  in  form  (20),  where  A\,  A\,  9,  and  9,  satisfy 
the  following  system: 

=^^^0*0*+-!-  sinOjj  ek>< , 

-H-^-(-^-cosO,  —  fjSinQjJeV  (/-I,  2).  (22) 

Here  we  use  the  designations 

%  -  //j#  sin  (8 — Pj)  —  Q,  cos  (0  —  pt), 

*1  *  ^ 4  +  I*i  c®»  (8  —  Pi)  —  «  sin  (8  —  p,))  Q,  —  //,  |/i,  sin  (8  —  -f- 

+  ocos(0  — pjlff. 
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%  —  HtRiinS  —  Q,cos6, 

X,  —  F,  +  AfQa  cos  6  —  sin  6. 


(23) 


that 


Let  us  now  deal  with  the  variables  z  and  z' ,  taking  into  account 

Al  -  AS*  ■  +  AA)  «M. 


Then  the  solution  for  system  (12)  is  written  as 

|  Z  —  A,COsQ,  +//jCO5(0— Px), 


|  z  —  A,  cos  Q,  cos  (0 
u  ■  AfcosQ,  +  cos  6, 


where 


"  » »  V  _ 

=  —  A,At  +  (i  |-ji-sin  Q,  +  cos  Q,) , 

'^~”'t('^COsQi  ~  ^sinQ*)  («-l.  2).  (25) 

Equations  (17)  and  (25)  are  reduced  to  the  standard  form.  We  will 
make  use  of  the  averaging  principle  for  their  solution  [2,  33.  After 
averaging  the  right  sides  of  equations  (17)  and  (25),  remaining  will  be 
the  terms  containing  only  the  variables  from  the  slow  arguments.  The 
equations  of  the  first  approximation  relative  to  the  variables  a>,  6,«,, 
Ax,  Vi.  At  and  have  the  following  form: 

«S  «  Q  =  const,  *  ( 26 ) 

J£-  »  _  A  A  —  H  g)  >l  COS  o  —  (a,  +  o)  sin  6], 

=  ~  H  l*iK  +  o)sin6  — tg  cos  6J  +  -^-j  +  . 


whe  re 


fl?  «»Vr(v?  — o*)*  +  4hio*. 


Taking  into  account  that  A  =  — — + —  and  A  «  o  ■*  jiA,  we  find 

A~J*  .  A~J »  (27' 

M-c-y,  ^  • 


7 


Substituting  (27)  in  the  first  equation  of  (26),  with  an  accuracy 
to  within  the  value  Of  the  order  of  p2 f  We  have 

*  *  **  2 <>•  —  [  *'«'*«?  8inP»H  q  sin26J.  (28) 

After  integrating  (28),  we  obtain 

®  -  ®,eip  |i  I - -M-U  + 


We  make  an  estimate  in  (19) 


^  |  sin  2&ft| . 


^-4  *3 

QA—C-J^rm*  ■  la  0, 

t 

We  note  that 


®«V 


A  — /i>0,  C—J't>0,  stop, 


■ >0. 
(Vj  —  o1)*  + 


It  is  obvious  that  in  order  for  <D  (/) -» 0,  it  is  sufficient  that  the  fol¬ 
lowing  condition  is  satisfied: 

(M-C-zjiWsinp,  ^  VC-J'J'  , 

,  ,1  «  ^  5  \  -J  -*■  / 

*«  fl  'I 

Since  (2A  —  C  —  J»f  —  (C  —  4)*  =  4  (A  —  C)  (A  —  /»),  the  following  relation  is 
satisfied  whenA>C: 

&A-C-Jfr>£-Jfr. 

Moreover,  if 

2/na*AlCT>4  >  mV*  I(v?  — a*)*  +  4A*o*], 

then  condition  (31)  will  be  satisfied. 

Thus,  condition  (3D  is  a  sufficient  condition  for  the  amplitude  of 
the  transverse  angular  motions  to  attenuate,  i.e.,  to  have  "true"  rota¬ 
tional  motion  around  the  X3  axis.  Weaker  sufficient  conditions  are  ob¬ 
tained  in  the  form 

I  A>C •  4  (32) 

\  2 maPhttrA  >  mV*  «v? -af  +  4A?o*l. 

For  the  case  A  =  0,  the  angular  speed  of  rotation  around  the  X3  axis  is 
constant  and  equals  Q.  As  can  be  seen  from  equations  (26),  in  the  case 
whe  re  the  oscillatory  motions  of  the  mass  at  the  first  damper 

will  be  attenuating  ones.  If 


8 


“ » + - - an^co, 

then  the  oscillatory  motions  will  not  be  attenuating  ones.  The  vibra¬ 
tional  motions  of  the  mass  at  the  second  damper  will  be  attenuating 
ones . 

When  A  *=  vlt  the  value  of  m,  is  not  constant  but  is  a  certain  func¬ 
tion  of  time.  In  this  case  the  solution  of  the  problem  gets  more  com¬ 
plicated.  Due  to  the  cumbersomeness  of  the  expressions  for  this  case, 
the  equations  of  the  first  approximation  are  not  presented. 

The  averaging  method  makes  it  possible  to  reduce  the  analysis  of 
angular  motions  of  the  body  supporting  the  rotor  and  damping  devices 
to  relatively  simple  equations  (26)  for  the  slowly  changing  amplitudes 
and  phases,  the  solution  of  which  can  be  easily  obtained  with  the  aid 
of  electronic  modeling  devices  or  digital  computers  and,  in  certain 
cases,  by  direct  integration.  The  analysis  of  motion  makes  it  possible 
to  point  out  its  following  special  features. 

The  movement  of  the  supporting  body  relative  to  the  center  of  masses 
will  be  close  to  the  Euler-Poinsot  motion,  at  which  the  angle  of  nuta¬ 
tion,  speed  of  proper  rotation,  and  the  amplitude  values  of  the  trans¬ 
verse  components  of  angular  velocities  of  the  supporting  body  change 
slowly  with  time.  The  energy  of  oscillations  of  elastically  suspended 
masses,  excited  by  the  revolution  of  the  supporting  body  and  rotor  with 
a  proper  selection  of  the  magnitude  of  masses  and  parameters  of  the  sus¬ 
pension  devices,  can  be  used  for  damping  nutation  oscillations  of  the 
system. 
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